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Abstract 

We derive the hydrodynamic equations for the supersolid and superhexatic 
phases of a neutral two-dimensional Bose fluid. We find, assuming that the 
normal part of the fluid is clamped to an underlying substrate, that both 
phases can sustain third-sound modes and that in the supersolid phase there 
are additional modes due to the superfluid motion of point defects (vacancies 
and interstitials). 

PACS numbers: 67.70. +n, 67.80.-s 
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I. INTRODUCTION 



Liquid helium ( 4 He) has a reputation for being the first substance in which one is able 
to observe many macroscopic quantum phenomena. In particular, it was the first system 
that could sustain superfluid flowjll and as a consequence display a number of amazing 
properties such as second sound, quantized vortices and the fountain effect. Furthermore, 
thin superfluid helium films were the first two-dimensional systems experimentally proven 
to undergo a Kosterlitz-Thouless transition to the normal state.! More recently it may have 
been observed!! that on weakly-binding substrates these films are the first-known spatially 
ordered superfluids.@ 

More precisely, measurements of the third- sound resonance frequency (which is propor- 
tional to the square root of the superfluid density) of submonolayer helium films on hy- 
drogen and deuterium substrates apparently indicate two independent Kosterlitz-Thouless 
transitions: the usual superfluid to normal transition at a temperature Tkt that obeys the 
expected universal jump relation,! and a second new transition at a temperature T c which 
is roughly 0.5 T^t for all coverages. The second transition appears as a sharp (but not 
discontinuous) rise or dip in the superfluid density depending on the substrate. 

In an attempt to explain these experimental results we have recently proposed that 
below the second critical temperature the superfluid helium film is in a spatially ordered 
phase exhibiting both off-diagonal (superfluid) and diagonal (hexatic) long-range order in the 
one-particle density matrix.! The main idea behind this proposal is that the hexatic to fluid 
transition is known to be a Kosterlitz-Thouless transition driven by disclination unbinding.@ 
(Disclinations are defects in the orientational order of a crystal created by the insertion 
or removal of a wedge of atoms, as shown in Fig. 1.) Therefore, our physical picture of 
the experiments is that at sufficiently low temperatures the film is in a superhexatic phase 
with only a dilute gas of bound vortices and bound disclinations present due to thermal 
fluctuations. For entropic reasons the disclinations then unbind at T c , leading to a transition 
from a superhexatic to a superfluid phase since the vortices remain bound at this transition 
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and the presence of free disclinations destroys the hexatic long-range order. At Tkt the 
vortices then also unbind and the film is finally forced into the normal liquid phase. 

Of course, to make sure that the above picture is qualitatively correct we must also 
consider the interaction between vortices and disclinations. This is even more pressing if 
one realizes that in a supersolid phase (where all disclination pairs are themselves bound 
into pairs or triples) this interaction is of long range and depends logarithmically on the 
distance between the two kinds of defects. Fortunately, it turns out that this is no longer 
true in the superhexatic phase due to the screening of the interaction by the surrounding 
gas of disclination pairs. A renormalization-group analysis actually shows that the vortex- 
disclination interaction is irrelevant and that the two separate Kosterlitz-Thouless transitions 
indeed survive. Nevertheless, the superfluid density is influenced in a non-universal way by 
the unbinding of the disclinations and Monte-Carlo simulations even show that on the basis 
of our hypothesis a rough qualitative agreement with the experiments of Chen and Mochel 
can be obtained.! 

However, to definitely identify the phase below T c more detailed information is needed. 
As a first step towards this goal we here present the two-dimensional hydrodynamic equa- 
tions of a superhexatic by describing the superhexatic as a supersolid with free dislocations 
(i.e. disclination pairs@). As a result of this approach we will also be able to consider the 
hydrodynamics of the supersolid phase, for which there is at present a renewed interest both 
in the context of Josephson-junction arrays! and solid 4 He.i Moreover, spatially ordered 
superfluid states have recently been proposed to be also relevant for the fractional quantum 
Hall effect ,0 since this effect can be understood as a condensation of composite bosons. § We 
therefore believe that the methods developed below might, if extended to bosons interacting 
with a Chern-Simons gauge field, also be used to obtain a description of the dynamics of 
such exotic quantum Hall states. 

We have organized the paper in the following manner. In Sec. [TI] we first consider the 
normal solid and hexatic phases by formulating a gauge theory that describes the phonons, 
the dislocations and the interaction between them. From this theory we then deduce in Sec. 
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1 1 J for both phases the dynamics of the appropriate hydro dynamic degrees of freedom. In 
Sec. [IV] we incorporate the effects of the additional superfluid order parameteiEl into the 
hydrodynamic equations derived in Sec. [TTI|and discuss the various long-wavelength modes in 
the supersolid and superhexatic phases obtained in this manner. We conclude in Sec. |VJ with 
a discussion on the possible relevance of our work to future experiments on submonolayer 
helium films and with a physical interpretation of our results. 



II. GAUGE THEORY OF PHONONS AND DISLOCATIONS 

In this section we will derive the long-wavelength (quantum) dynamics of the solid and 
hexatic phases. The discussion closely follows work by Kleinert,^ save that we will not 
include higher gradient elasticity. This leads to a considerable simplification of the theory 
but implies that we cannot properly treat the dynamics of the disclinations. Fortunately, for 
our purposes only the dynamics of the dislocations is of importance and this simplification 
is justified. 



A. Solid 



In the case of an isotropic crystal, the action for the displacement field Ui(x, r) in the 
presence of a pair of dislocations is given byE! 



S[ui] = j o drj dx i^(d T Ui - A) 2 + fi [uij - ^L_n± j + _ p u y j, . ( i ) 

where u%j = (diUj + djUi)/2 is the strain tensor, p and A are the usual Lame coefficientsQ 
and p is the average mass density. The unphysical (and singular) contributions arising from 
the multivaluedness of Ui(x, r) are compensated by the quantities f3j and (also known as 
the 'plastic distortion'). Their relationship to the defects is best explained by the Volterra 
construction^ Let C be a small loop bounding a section of two dimensional crystal that is 
excised from the whole (cf. Fig. 2). The edges of the loop are drawn together and form a 
line C. This line may be time dependent, and its definition is not unique. However, the 



topological defects (i.e. two dislocations with opposite Burgers' vectors) associated with the 
distortion of the surface are always located at the endpoints of £. If ±B are the Burgers' 
vectors of the dislocations constituting the pair and if v is their velocity then = Si(£)Bj 
and Pj = —Vi5i(C)Bj. The delta function 5{(C) is singular on the time-dependent Volterra 
cutting line £ of the dislocations and is directed along the normal vector. If the cutting line 
£ is parameterized by x(s, r) with < s < 1, this means mathematically that 

8i(C) = -eyfds 9Xji d S s ,T) 5(x-x(s,t)) , (2) 

where is the two-dimensional antisymmetric tensor. Note that the dislocations are as- 
sumed to be able to move freely, without any friction, through the crystal because the 
equations of motion for the displacement field allow for time-dependent solutions that pre- 
cisely correspond to such evolutions of the crystal.B We will come back to the issue of friction 
in Sec. |T| when we consider the effects of dissipation. 

We now first perform a Hubbard-Stratonovich transformation by introducing the auxil- 
iary variable p (representing the momentum density) and adding the quadratic term 



rhP r 1 

dr J dx — (p* - ip(d T Ui - Pi))' 



2p 

to the action, which may now be rewritten as 



S[pi,Ui] = J dr J dx + fi (uij - 



2p r V 2 J 

+ ^(uu - Pa) 2 - ipi(d T Ui - Pi)\ . (3) 

Integrating out p would return the original action up to an unimportant constant. In a 
similar manner we then also introduce the symmetric stress tensor er^-, to decouple the 
terms quadratic in the strain. This results in 



Sfo^M = fdrjdx [ | + -L (4 - ^-al) 

Pij Pj 



2p 4/x 

ipi{d T Ui - pi) + iaij [ Uij - ^ - ^' Jl \ \ , (4) 
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with v = A/(2/i + A). The partition function is now given by the functional integral 

Z = J d[pi] J d[aij\ J d[ui\ exp j-^^, cr^, ujj , (5) 

where the integration over is only over the symmetrical part since we have not included 
higher gradient elasticity. 

We can now perform the integration over the displacement field. Because the action is 
linear in Ui this simply leads to the constraint 

d T pj = didij . (6) 

This constraint can be automatically satisfied if we introduce the vector field A,- and the 
tensor field by setting = t ik d k A,j + t ki d T A k j and pj = e k idiA k j. Substituting these 
relations into the action we find that the interaction between the gauge fields (i.e. the 
phonons) and the dislocations is given by 

rh/3 r 

S int [A i:j , Aj] = J dr J dx {-iA^ + iA ij J ij } , (7) 

where after several partial integrations the unphysical singularities of $ and have disap- 
peared and only the dislocation density and the dislocation current density remain. Intro- 
ducing also the function 5(V) = S(x(l, r)) — 5(x(0, r)), which denotes the difference between 
a delta function at one end of the cutting line C and a delta function at the other end, these 
densities and currents can conveniently be written as otj = 5(V)Bj and J^- = —Vi5(V)Bj, 
respectively. As a direct consequence of the above definitions they obey the conservation 
law 

d T ocj = diJij . (8) 

In addition, the dynamics of the phonons is determined by the remaining quadratic terms 
in the action which expressed in terms of the gauge fields Aj and Aij yield 



with a-jj equal to €i k (d k Aj — d T A k j). Comparing this result with Eq. flU) we observe that 
the stress and the physical part of the strain u^ hys = Uij — (fy + (3ji)/2 are related by 
crjj = 2/j l Utf hys + X5ijUkk Va and therefore by 



We will have need of the latter relation in Sec. |ITT| , when we discuss hydrodynamics. A more 
formal way to justify it is to add to the action S[ui] a source term 



and perform the same manipulations as before. We then find that the source Kij indeed 
couples linearly to the right-hand side of Eq. ([U]) . 

Following Kleinert, we now notice that the above theory has a gauge symmetry as a result 
of the fact that the gauge fields Ai and Aij are not uniquely determined if the stresses exp- 
and momenta Pi are known. Indeed, <7y and Pi are invariant under the gauge transformation 
Ai —>■ Ai + <9 T Aj and A^ —>■ A^ + diAj. Hence So[A{j, Aj] is also invariant. Moreover, due to 
the conservation law in Eq. (|8]), the interaction Si n t[Aij, Aj] is invariant too. 

To calculate the partition function we therefore need some gauge-fixing procedure. The 
symmetry of requires that 



We would now like to write the gauge fields as the appropriate derivatives of unconstrained 
fields. Using the above gauge symmetry we can always take Ai = eijdjX and An = 0. This, 
however, does not completely fix the gauge because these conditions are still invariant under 
the smaller group of transformations \ — * % + 9 T A and Ay — > Ay + (^(e^c^A). To see more 
clearly the consequences of this residual symmetry we expand A^ into its longitudinal and 
transverse components (with respect to both indices), i.e. 




(10) 




^ij&ij — 9jAj — d T [Ajj) — . 



(11) 



Aij = di(djA ) + di(e jk d k A LT ) + e lk d k (djA TL ) + e lfc <9 fc (e j7 «9^ TT ) , 



(12) 
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where we have introduced four new fields. The tracelessness of A%j can then be fulfilled 
by taking A LL = —A TT . In addition, the residual gauge symmetry can now be written as 
X —>■ X + 9 T A and A LT — ► A LT + A. This shows that instead of the fields x an d A LT we must 
use the gauge-invariant field x' = X ~ d T A LT together with A as integration variables. The 
associated change of measure can be incorporated in the normalization and the same is true 
for the 'volume' / d[A] of the residual gauge group because the action is gauge invariant and 
therefore cannot depend on A. After this gauge-fixing procedure the partition function thus 
becomes 

Z = J d[A TT ] J d[A TL ] J d[ X '\ exp {-i (S [A TT , A TL , X ') + S mt [A TT , A TL , X '))] . (13) 

Note that we are left with three physical degrees of freedom, which is the correct number in 
two dimensions since (Ty and p, contain in principle a total of five degrees of freedom but we 
have two constraints in Eq. (^]). Note also that the transformation from and pi to A TT , 
A TL and x' is a linear one so that the Jacobian involved in the calculation of the partition 
function is simply an unimportant constant. In particular, the stress is given by 

= tiktjtdkdix' + d T (didjA TL + e ik d k d j A TT + e jk d k diA TT ^ , (14) 

which is manifestly symmetric in % and j. 

A straightforward calculation now shows that the free part of the action is 



S [A TT ,A TL , X '} = /""dT J dx {^(drd'A^r + j-id^A 



2 aTT\2 



- T,Th {a ^ AT ^} ■ (15) 

It contains four modes: The part involving A TT has a pair of modes (corresponding to ±fc) 
with uj 2 = fik 2 /p. These modes therefore represent the transverse phonons with a speed of 



sound yyu/p. The part involving A TL and x' nas another pair of modes with a dispersion 
obeying uj 2 = (2/x + X)k 2 /p. These represent the longitudinal phonons with a speed of 
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sound y (2/i + X)/p. Interestingly, these results can be understood much more easily if we 
introduce the field 

X " = X ' ~ vdrA TL , (16) 
since then the above action becomes 

+ ^ m + ^A^f + j^—^f } , (17) 

so that the fields are completely uncoupled. Notice that the x" field has no kinetic term, 
which explains why we obtained above only four modes instead of six, as might have been 
expected in first instance. 

— # 

Furthermore, if we introduce the usual Burgers' field b(x } r) for the total dislocation 
density, which is nothing more than the sum of the density over all dislocation pairs, then 
the interaction with the dislocations acquires the form 

S int [A TT , A TL , x'} = drjdxt {\'<, /),!>, - A TT d T dA} , (18) 

where we have made use of Eq. (|j) to express the longitudinal part of the current density 
Jij in terms of the time derivative of hi. Decomposing b into its transverse and longitudinal 
parts, i.e. hi = d{b L + eijdjb T , the interaction finally becomes 

Sint \A TT , A TL , x') = drjdxi { X 'd 2 b T - A TT d T d 2 b L } . (19) 

The total action S = So + Si nt reduces for time-independent configurations to the one we 
previously used for a discussion of the critical properties of the superhexatici Integrating 
out the fields A TT and x' we can now fi n d the time- dependent interaction among the dis- 
locations. Physically, these interactions are thus associated with phonon exchange and the 
time dependence arises due to the finite speeds of sound. This picture also explains why the 
effective action for x' contains just one pair of modes: The self- interaction of the transverse 
dislocation density can only be mediated by longitudinal phonons. 
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B. Hexatic 



Up to this point the dislocation density has not been an independent dynamical variable, 
since we have specified the positions of the dislocations at all times and thus neglected the 
influence of the phonon dynamics on their motion. However, to describe the hexatic phase 
we want to integrate out the dislocations in the plasma (or continuous) approximation.@ 
For that we need the free action of the field b. Here we can again make use of the results 
obtained by Kleinert, who showed that the energy associated with the nonlinear stresses at 
the heart of the defect can be lumped into a 'core contribution' to the action.lll0 In our 
notation this contribution becomes 




This action represents free propagation of the dislocation density fluctuations which, as 
mentioned previously, are permitted by the classical equations of motionEl and neglects 
dissipative coupling of the dislocation cores to the phonons. 

Integrating out the Burgers' field using Eqs. (|TJ) and (p0|), we obtain the following results. 
The effective action for A TT becomes 

+ 2k< a ^ TT )(0 +1 )" ( ^'' 4lT) }- (21) 

As shown in Fig. 3a, it contains two pairs of modes which for k 2 ^> p/2E c all have a dispersion 
obeying uj 2 ~ pk 2 / p. However, for k 2 <C p/2E c , one pair of modes has a dispersion uj 2 ~ 
p 2 /E c p + 2pk 2 /p with a gap whereas the other pair of modes is gapless with uj 2 ~ 2E c k 4: /p. 
This is consistent with our expectation that in the hexatic phase there should only be one 
pair of transverse gapless modes with a softer dispersion than that of the transverse phonon 
modes in a true solid. 
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The effective action for \' an d A TL in the hexatic phase is 

^W^ 1 )" (22) 



Integrating out also A TL we finally arrive at the effective action for \' ■ It reads 



-i 



2 / o9 o 2 \ -1 



(<WV)> (23) 

and also contains two pairs of modes (cf. Fig. 3b). For k 2 3> p/2E c we recover of course the 
ordinary sound dispersions uo 2 ~ 2fik 2 /(p(l — v)) = (2/i + X)k 2 /p and c<j 2 ~ pk 2 / p. However, 
in the limit k 2 <C p/lEc these evolve into a pair of gapped modes with u 2 ~ 2p 2 / (E c p(l — u)) 
and a pair of propagating modes with oj 2 ~ 2/x(l + v)k 2 /p, respectively. Clearly, the same 
mode structure is also present in the effective action for A TL (obtained by integrating out 
X 1 instead of A TL ) which indicates that in the hexatic phase the longitudinal velocity is 



renormalized downwards to y2p(l + u)/p. 

III. HYDRODYNAMICS OF SPATIALLY ORDERED PHASES 

We now turn to the linear hydrodynamics of the solid and hexatic phases that follows 
from the theory presented above. For the sake of clarity, and because it will turn out to be 
less important for our purposes, we will not discuss temperature fluctuations in the following. 
However, having derived the relevant energy densities in Sees. |11 A| and [11 B| it is in principle 
straightforward to include temperature fluctuations in our theory and, in particular, to arrive 
at the extension of the hydrodynamic equations presented below that is required if one wants 
to consider also the hydrodynamic mode due to energy conservation. After the equations 
of motion for the hydrodynamic variables are determined, we can find the propagating and 
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diffusive modes. This is done as before, by Fourier transforming the equations of motion 
and determining the dispersion u(k). Propagating modes appear as complex roots of a 
characteristic equation and will always occur in pairs. Each physically distinct propagating 
excitation such as longitudinal or transverse sound corresponds therefore to two roots or 
modes. 

We start by considering the mass-density fluctuation Sp above the average mass density 
p and initially neglect the possible presence of vacancies and interstitials. To lowest order 
in the strain, the density fluctuation equals —pu^ hys so up to that order we obtain 

d T 6p = -pd T ui hys = -P—fan , (24) 

if we make use of Eq. ( |TH| ) to relate the stress and the strain. Using also the decomposition 
Cii = d 2 x' + d T d 2 A TL from Eq. (|T4P we may write this as a pair of continuity equations 

d T 5p = d i9 l , (25a) 

d T gf = durfj , (25b) 
where the longitudinal momentum density is given by 

9 '=-^TA) S ^' + ^ Ti) (26) 

and the diagonal part of the stress tensor by 

rf = -^Ah-M + d rA TL ). (27) 



In the absence of defects, the hydrodynamic quantities gf and ir? are precisely equal to 
the longitudinal part of Pi and the diagonal part of a^, respectively. This can be seen in 
the following manner. In an ideal solid we have no defects, and variation of the action in 
Eq. fllTD gives x" = or x' = ud T A TL , which we may use to eliminate \' ■ The equation of 
motion for A TL generated in this way is 



>2 atl _ 2/i + A 2 TL 



d 2 A iL = _ zn _ llQ 2 A rL (2g) 
12 



If we substitute this back into the definitions of gf and irfj we obtain the longitudinal part 
of Pi and the diagonal part of as given in section ||. In the presence of defects with their 
own dynamics this is no longer true, since the dislocation density couples to the gauge fields 
and alters the equations of motion. To avoid confusion about this point we have, therefore, 
introduced a new notation for the hydrodynamic momentum density and stress tensor which 
we will use for the rest of the paper. 

We also note that the above equations are not Galilean invariant and are therefore only 
valid in a specific reference frame. This is a result of the fact that the gauge theory of Sec. [II 
has implicitly used the existence of an ideal lattice with respect to which the displacements 
u(x, t) are defined. El Hence, the prefered reference frame corresponds to that frame in which 
this ideal lattice is at rest. This is the case for all the hydrodynamic equations presented 
below. 

In this ideal solid without interstitials or vacancies the pressure fluctuation (following 
from ttP- = —5ij5p) equals 



and the mass-density fluctuation becomes 



Substituting the equation of motion Eq. (|28| ) for A TL into Eq. (|29|) , we obtain the desired 
constitutive equation 

5p = 2 -^6p . (31) 

Together with the hydrodynamic equations (p5|) this correctly leads to the sound equation 

d *S p = _lE±±& S p = -^& S p y (32) 

with C|| the longitudinal sound velocity. 

However, as stressed first by Martin, Parodi, and Pershanlil and again by Zippelius, 
Halperin, and NelsonS we are not in general allowed to assume that the crystal is ideal, 
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without vacancies or interstitials. We must include the effects of (long-wavelength) fluctu- 
ations in the net defect density n A , which is defined as the density of vacancies minus the 
density of interstitials. To do so we can make use of the fact that in a hexagonal system these 
defects can be regarded as a 'bound state' of three dislocations with radial Burgers' vectors 
pointing symmetrically outward (interstitial) or inward (vacancy) @ This is illustrated for 
an interstitial in Fig. 4. As a result the interaction of the net defect density with the gauge 
fields is given by 

rh/3 

ir I dx % 

\i 

where Vq denotes the area deficit induced by a defect in an otherwise perfect crystal and 
7a = ^c^V /2c 2 ± . We can verify this result by noting that in the static case (and n A — > in A 
because of our conventions in the imaginary time formalism of Sec. [TT]) the Euler-Lagrange 
equation for the Airy stress function, following from the action in Eq. flT5] ) together with 
the above interaction, becomes d 2 x = 2(/x + A)VonA which correctly leads to / dx u^ hys = 
V J dx n A . Furthermore, the free action of the defects becomes (cf. Eq. (|20|)) 



S mt [A TT , A TL , x'} = r dr ( dx t^n A d 2 X ' , (33) 
jo J u 



rhp r E A fpd 2 \ 
S [n A } = drj dx —^-n A I + 1 I n A , (34) 

with E A of order E c Vq. 

Redoing our calculations with n A non-zero, we find that n A displaces the \' field. There- 
fore x" i n Eq. ( |TED is also non-zero. Moreover, we now obtain instead of Eq. ( |3"2| ) the coupled 
set of equations 



27a \ d 2 nA + i 7aA 
E A ji J E A p 2 



d 2 n A = - c {[l + -^\ d 2 n A + i-^d 2 5p , (35b) 



for the longitudinal degrees of freedom. Note that the density fluctuation 5p receives a 
contribution from both the lattice vibrations as well as from the net defect density, since 

S P = -\ dr d 2 A TL + ^n A . (36) 
14 



As a result the longitudinal momentum density has also two contributions 

g L = -\d id *A TL - 2 -^Jt , (37) 

c ll c ll 

where J L is the longitudinal part of the net defect current density obeying the continuity 
equation d T n& = —diJf. 

This almost completes our discussion of the hydrodynamical description (without dissi- 
pation) of the solid phase. However, we have not yet obtained the transverse modes. From 
our expressions for the strain tensor Uij one can easily show that in the solid phase the 
transverse part of the displacement field is given by 

uj = ^-e ij d j (d T A TT ) = e ijdj(d T A TT ) , (38) 

where c± is the transverse speed of sound. Hence, the transverse dynamics of the lattice 
is solely determined by the transverse phonons and we have the additional hydrodynamic 
equation 

d 2 T A TT = -cld 2 A TT , (39) 

which is completely uncoupled from the previous ones and in particular does not depend on 
the net defect density ha- Moreover, if we introduce the standard hexatic order parameter 
#6, which is equal to the local bond angle and may therefore be written as 

tf 6 = \e i3 d lU] = -—^ T (d T d 2 A TT ) , (40) 

the equation for A TT is equivalent to 

%A = -c 2 ± d 2 V 6 , (41) 

so that i? 6 can also be used to describe the transverse phonons. 

From Eq. (j38|) we also find that the transverse part of the momentum density is given 

by 

91 = ^(0^). (42) 

c ± 
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Therefore the stress tensor has the nondiagonal contribution 

1 



7T 



r = ~^m TT ) = 2pcle ij A , (43) 



and both the longitudinal as well as the transverse hydrodynamic equations in the solid 
phase can be summarized by 

d5p 



Of 



-V • g , (44a) 



-JL = - c 2 V5p - 7 A Vn A + 2pc 2 ± V x tf 6 , (44b) 



d °° '-Vxj, (44c) 



dt 2p 
dn A 



V ■ J , (44d) 



dt 
df 

— = -c A Vn A + iV5p , (44e) 

after a transformation to real time, which in this case not only means that r — > it but also 
g — > ig and n A — > in&. Moreover, note that the constants c, 7a, ca, and 7 should here 
be interpreted as renormalized quantities which are determined in terms of the microscopic 
parameters of our gauge theory by a comparison with Eq. (pop. 

Now we are ready to discuss dissipation. In principle dissipation has already been in- 
cluded because there is a coupling between the net defect density tta and the phonons. 
Hence if for example an interstitial were, in a discrete picture, to tunnel from one loca- 
tion to another there would be a 'shake up' of the phonon field. However, if we treat tta 
as a smooth continuously varying field, the action in Eq. ( |34"|) is quadratic and the bilin- 
ear coupling n&d 2 x' produces only mixing of the collective modes but no real dissipation. 
Therefore, we choose to include effective dissipation in the same manner as explained in 
detail by Zippelius, Halperin and Nelson.0 Using their notation we first of all add to the 
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right-hand side of Eq. (|44b| ) the terms {rjS/ 2 g + CV(V ■ g))/p associated with the dissipative 
part of the stress tensor 7Tjj and representing viscous diffusion of the momentum density. 

Next the question arises how we need to modify Eq. (|44e|) . This equation is a result of the 
fact that we have allowed the dislocations, and hence the interstitials and vacancies, to move 
freely through the lattice and used Eq. fl3H|) for the free action of the defects. If the defects 
effectively experience friction (for example due to their interaction with the phonons), then 
it is more appropriate to add a Leggett friction term0 and use 

S [n A ] = j Q drj dx -^n A ( ^ + i^d T + lj n A (45) 

instead. The dispersions then indeed obey u)^ ~ ±cj_k — iC,k 2 /2 at long wavelengths, and 
we must add the term £V(V ■ J) to the right-hand side of Eq. (|44e ). If we further assume 



that the transverse part of the defect current density behaves as in a gas and simply diffuses 
to zero with a diffusion constant k, we obtain in total 

85p 



dt 



-V • g , (46a) 



?2 = --V5p - 7 AVn A + 2pciV x ^ 6 + ^V 2 g + -V(V • g) , (46b) 
ot p p p 

? = -Vxg, (46c) 

dt 2p y v ; 

<£ = -V./, (46d) 

-4Vn A + 7V5p + kV 2 J + «eV(V • J) , (46e) 



dJ_ 
~dt 



with B = pdp/dp\ nAt T = pc 2 the appropriate isothermal bulk modulus in view of the fact 
that the pressure is a function of both the particle density as well as the net defect density. 
From thermodynamics we therefore also conclude that 7a = dp/dn A \ Pt T- 

It is interesting to point out that these hydrodynamic equations differ from the results 
obtained by Zippelius, Halperin, and Nelson. In particular, their Eq. (3.32) differs from our 
Eq. ( [46e|) and reads 
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J=-T A V\^-- lA 6pj . (47) 

The difference can easily be traced back to the fact that Zippelius, Halperin, and Nelson 
assume on phenomenological grounds that the dynamics of the net defect density is purely 
diffusive. Indeed, we exactly reproduce their results if we use in our calculation a free action 
of the form 



S [n A ] = dr J dx ( — + lj n A (48) 

instead of Eq. (|34]) . We can therefore consider the hydrodynamic equations of Zippelius, 



Halperin and Nelson as the overdamped (or classical) limit of our Eq. (f46|). Clearly, Klein- 
ert's more microscopic approach does not lead to purely diffusive but in first instance to 
propagating behavior of the defects, which is appropriate for the quantum crystals of in- 
terest in Sec. [TV[ We now turn to the modification of the above results in the hexatic 
phase. 

In the hexatic phase there are free dislocations present and x" replaces n A as the appro- 
priate dynamical degree of freedom. To see most clearly how this comes about we will work 
perturbatively in 1/E C . Up to first order in 1/ E c the effective action for x" is 



which upon Fourier transformation displays two modes with uo 2 = c\k 2 + 2// 2 (l + v)j[E c p). 
So in the limit E c — > oo (which physically means that we are looking at the nonhydrodynamic 
regime k 2 ^> ji/2E c ) we approximately have 

d 2 rX " = -c\d 2 x" , (49) 
whereas the equations of motion for A TL and A TT are 

d 2 A TL = -c\V 2 A TL + Y^drX" (50) 

and 



d T A = —c ± V A , (51) 



respectively. For the mass-density fluctuation we now find 



f ' = H^ lTi -2o£*j 8Y (52) 



and for the stress tensor 



4 - j^" + 5<^V] - h- ^Y'} • (53) 

Putting all this together we obtain in first instance the following set of hydrodynamic 
equations for the hexatic phase 

d T 5p = V • g L , (54a) 



d T g L = -V5p = -c\V5p - ^—^(W ) , (54b) 



dh" = -clvV , (54c) 

d 2 T A TT = -clV 2 A TT . (54d) 
Combining the longitudinal and transverse parts as before, this equals 

d T 5p — V • g , (55a) 



d T g = -c\V5p - 2(1 1 +z/) V(VV) + 2pciV x tf 6 , (55b) 



dh" = ~cW , (55c) 



d T A = ~^-Vxg, (55d) 
2p 

and clearly reduces to the hydrodynamic equations for the ideal crystal if we put x" — 0- 
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We now have to consider how the above picture changes for a finite value of E c . Here 
we can use the results of Sec. [11 BL In the hydrodynamic regime k 2 <C p/2E c we saw that 



the transverse speed of sound was renormalized to zero, because we found the quadratic 
(particle-like) dispersion u 2 = 2E c k 4 / p. As a result we have for the transverse part of the 
hydrodynamic equations 

d 2 A TT = ^ d A A TT ) (56) 



which implies that in the right-hand side of Eq. (|55b|) we must replace 2pc 2 L 'V x $6 by 
-AE c e z x V(V 2 $ 6 ). This gives 

d T f = -AE c e z x V(V 2 tf 6 ) , (57a) 

d T A = -^V x f , (57b) 

which is in complete agreement with Zippelius, Halperin, and Nelson if we identify the Frank 
constant with 8E C . 

For the longitudinal part we need to analyze the dynamics of x" an d A TL . A straight- 
forward calculation shows that the effective action for these fields contains precisely the 
four modes already found in Sec. [11 H. The propagating modes with uj 2 = 2/i(l + v)k 2 j p 



obey x' = x" + vd T A TL = and are therefore indeed associated with density fluctuations 
proportional to d 2 x" ■ We thus need to use a renormalized longitudinal speed of sound equal 
to 



c= Ml + ») = PP2P + 2X 

V P V P 2^ + A y 1 

that is always smaller than the longitudinal speed of sound in the solid phase. In fact, 

this actually exhausts the longitudinal hydrodynamic modes since the other modes in the 

effective action for x" an d A TL are gapped. As a result we now obtain in real time the 

following set of hydrodynamic equations for the hexatic phase 

IT = < 59a > 
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^ = -c 2 VSp - ^ S, x V(V 2 tf 6 ) , (59b) 
f = ^Vx ? , (59c) 

not including dissipation. 

To include dissipation we again follow Zippelius, Halperin and Nelson and add to the 
right-hand side of Eq. ( [59 bp the terms {rjS/ 2 g + (V(V ■ g))/p- However, we do not add 
the term nW 2q d§ to the right-hand side of Eq. ([59c]) because, just as in the solid phase, the 
dissipation of the transverse modes is already accounted for in the term r]S7 2 g that is added 
to the momentum equation. Put differently, a term of the form kV 2, &q can be absorbed by 
an appropriate redefinition of Ka and r\. Again introducing the isothermal bulk modulus 
B = pdp/dp\ T = pc 2 we then find 

d5p 



Of 



-V • g , (60a) 



% = --VSp - ^e z x V(V 2 tf 6 ) + V -V 2 g + ^V(V • g) , (60b) 
at p I p p 

— = -V X g (60c) 

as our final result for the hexatic phase. Apart from the absence of a dissipative term in Eq. 
(|60cj) it agrees with the findings of Zippelius, Halperin and Nelson and therefore contains the 
same mode structure as derived in that paper. For completeness sake, we mention however 
that the equations of motion for the hexatic order parameter $6 can be derived from an 
effective action 

r n/3 r I / d 2 \ 
S e "\M = j o drj dx -A Up-£ + 4t V d T - K A d 2 \ A , (61) 

that can easily be understood physically: The first term on the right-hand side corresponds 
to the kinetic energy / dxp(d T u) 2 /2 of the displacement field. The second term is a Leggett 
friction term and the last term corresponds to the usual Frank energy, which is responsible 
for the fact that the hexatic to liquid transition is of the Kosterlitz-Thouless type. 
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IV. HYDRODYNAMICS OF SUPERFLUID PHASES 



Having arrived at the hydro dynamic equations for the solid and hexatic phases, our next 
objective is to incorporate the effects of the additional hydrodynamic degree of freedom 
associated with the phase of the superfluid order parameter. Fortunately, from the micro- 
scopic theories developed for superfluid liquidsil and gasesil it is well known how we should 
proceed to obtain the hydrodynamic (two-fluid) equations for the superfluid phases starting 
from the equations for the normal phase. The procedure consists in principle of four steps. 
First, the total (average) density p of the system is split up into a normal density p n and a 
superfluid density p s . In general these densities are tensors of second rank, but for systems 
with hexagonal symmetry which are of interest here they are proportional to the identity 
5ij and can be considered as scalars. Second, the total momentum density g is similarly 
split up into a normal component p n v n and a superfluid component p s v s with a superfluid 
velocity that is purely longitudinal (V x v s = 0). Third, for an effectively isotropic system 
the dissipative terms in the momentum equation must be generalized to 

r]V 2 v n + Ci— V(V • {v s - v n )) + C 2 V(V • v n ) . 
P 

Finally, we must add the dynamics of the superfluid velocity, which is basically determined 
from the Josephson relation and reads 

^§ = -^VSp + Cs-V(V • (v s - v n )) + C 4 V(V • v n ) , (62) 
at p 2 p 

where B = p 2 dpjdp\T is again the isothermal bulk modulus and p is the chemical potential 
per unit mass. We again leave out temperature fluctuations since we are primarily interested 
in third-sound modes, for which these fluctuations are (at least qualitatively) unimportant. 



A. Supersolid 

To apply the above procedure to Eq. ( pEBP we must realize that we are here in fact 
already dealing with a two-fluid hydrodynamics. We must therefore not only split up the 
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total momentum density g into a normal and a superfluid component but also the net defect 
current, i.e. J = J n + J s . Moreover, we have to account for the fact that the chemical 
potential, just like the pressure, is a function of the particle density and the net defect 
density. In this manner we arrive at the following hydrodynamic equations 

d5p 



dt 

B. 



-V • g , (63a) 



t^- = V£p - 7a Vn A + 2pciV x $ 6 

at p 

+ r]W 2 v n + Ci-V(V • (v s - v n )) + C 2 V(V • v n ) , (63b) 
P 

(63c) 

^ = -^V5p + [3 A Vn A + (AV(V ■ (v s - v n )) + C 4 V(V • v n ) , (63d) 

-^ = -V-J, (68b) 
df -> 

— = -4Vn A + 7 V5p + kV 2 J n + £iV(V • J s ) + £ 2 V(V • J n ) , (63f) 

^ = -^Vn A + /?p s V5p + £ 3 -V(V ■ J s ) + £ 4 V(V • J n ) , (63g) 
or p z p 

with /?a = —dfi/dn A \ Pj T- These represent nine equations for the nine unknown functions 
Sp, v n , v s , i? 6 , n A , J n and J s . 

Although a complete analysis of the various hydrodynamic modes is now possible, we 
will consider here only the situation which is most relevant to experiments, namely that the 
normal part of the two-dimensional system is clamped to an underlying substrate. As a 
result we have v n = J n = and Eqs. ( |63b| ) and ( |631| ) determining the normal properties of 
the supersolid are no longer valid. The hydrodynamic equations therefore reduce to 

^ = ^V 2 5p - p APs Vn A + Cs^|(V 2 5p) , (64a) 
at p 2 p at 
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d n A B A p s 2 2 p s d 2 

-£T = — 2~V n A - (3p s V 5p + £ 3 — — (V n A ) . (64b) 
ot p z p at 

They contain two pairs of propagating modes, which in the limit of a small coupling constant 
(3 <C BB A / (3 A p i essentially correspond to a pair of third-sound modes with 5p unequal to 
zero but a constant net defect density and a pair of modes with an oscillating net defect 
density. 

One might have expected that the coupling of a superffuid density to a propagating 
defect density would have resulted in one pair of gapped excitations and one pair of gapless 
excitations instead. Consider, for example, two identical superfluid layers. If the layers are 
uncoupled the dynamics of the phases *&i and $2 of the layers is determined by the action 



^layers 



[A, A 



fdr J dx [^(d T A) 2 + I (WO 2 + ^{d T A? + f (W 2 ) 2 } , (65) 



which clearly has two pairs of gapless (third-sound) modes, one pair for each superfluid. If 
we couple the order parameters by allowing the particles to tunnel with an amplitude —J/p 
from one layer to the other we must add a Josephson coupling 



S-, 



tunnel 



[tfi, tf 2 ] = - fdr f dx Jcos(tfi - # 2 ) (66) 



to this action. The hydrodynamics modes couple to form two in-phase and two out-of-phase 
excitations. The modes with d\ and i?2 oscillating out of phase get gapped (i.e. uj 2 ~ BJ '/ p 2 
for k 2 J I ' p s ) and only the modes with d\ and i?2 oscillating in phase remain gapless. Yet 
in Eq. ( |64] ) we find only gapless modes. 

This paradox can be resolved by noting that we have made the standard assumptionlil 
that both the total number of particles and the net number of defects is conserved. Hence, 
after an atom has tunneled from a lattice site to the position of a vacancy, a new vacancy 
is created near the original site of the atom. The analogous process for the two coupled 
superfluid layers in not simply tunneling of individual atoms from one layer to another, but 
rather the exchange of a pair of atoms in different layers, returning the system to its original 
state. Such a process is not a Josephson coupling and therefore the modes remain gapless. 
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The existence of separate conservation laws for the particle and defect density thus allows 
in principle two separate broken symmetries. 

We also note in passing that the third-sound modes in Eq. ( p^ ) are not present in 
the hydrodynamic equations proposed by Andreev and Lifshitzil and considered in more 
detail by Liu.@ This is a result of the fact that these authors use a somewhat different 
physical picture for the supersolid phase: They assume that the superffuid current density 
is carried by (Bose condensed) defects and that the normal current density is solely due to 
lattice vibrations. Hence if we take v n = 0, which in their context means that du/dt = 0, 
only transport of defects is possible and only the latter two modes survive. However, as a 
consequence of their picture the hydrodynamic equations in the (normal) solid phase describe 
only longitudinal and transverse sound modes in an ideal lattice and do not include the effect 
of vacancies or interstitials. As explained above this is incorrect in principle and one should 
at least also allow for a normal current density due to the motion of defects. In addition, we 
have seen in Sec. |T| that even in the presence of defects the density fluctuations are equal 
to —pu^ hys . It is therefore perfectly reasonable that if there is superffuid mass transport 
possible in the solid, it can be caused both by the motion of defects and by lattice vibrations. 
Indeed, as an existence proof of this latter possibility we can for instance consider superffuid 
4 He in a weak periodic and commensurate potential, which is clearly a supersolid without 
defects. 

While it is generically possible to have both density and defect superffuid modes, we 
might expect however, for realistic films on realistic substrates, that in a supersolid it may be 
harder for particles to perform ring exchanges0'il than for vacancies to exchange positions. 
Thus, a priori, we might expect the effective superffuid stiffness for the density fluctuations 
to be smaller than that of the vacancies, perhaps to the point where the former is entirely 
absent. 
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B. Superhexatic 



We next turn to the superhexatic phase. In a similar manner as in Sec. IV A we obtain 
from Eq. (|5D|) the full set of hydrodynamic equations 

d5p 



dt 



-V • g , (67a) 



f = --V5p-^e 2 xV(V¥ 6 
at p 2 

P 



+ rfJ 2 v n + (i-V(V ■ (v s - v n )) + ( 2 V(V • v n ) , (67b) 



^ = 2^' (67C) 

^ = V5p + Ca- V(V ■ fa - if n )) + C 4 V(V ■ v n ) , (67d) 

that leads to the usual two-fluid hydrodynamics of a superfluid if we omit Eq. (|67c[) and 
put i9 6 = 0. Therefore these equations allow for first and second sound,il and for a pair of 
transverse modes involving and $6 which are either dispersive or propagating depending 
on the sign of A = K^/Ap — (rj/p n ) 2 : If A < we have two purely dispersive modes with 
c<j ± = —i(i]/p n ± V - A)Ar/2, wheras if A > we have two propagating modes and the 
particle-like dispersion lu^ 1 = iv^A k 2 /2 — i(rj/p n )k 2 /2. However, considering again the case 
v n = the hydrodynamic equations now simply reduce to 



which contains only a pair of third-sound modes with the velocity c 3 = ^Bp s /p 2 and the 
diffusion constant D 3 = (,zp s / p. 



V. CONCLUSIONS AND DISCUSSION 



In this paper we have derived the hydrodynamic equations for the supersolid and super- 
hexatic phases of a neutral two-dimensional Bose fluid. For the supersolid these equations 
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are rather complex, since they incorporate the effects of defect motion and lattice vibrations 
on both the normal and superfluid parts of the momentum density. Our physical picture 
for the influence on the superfluid part is roughly speaking that in a mean-field theory the 
condensate wavefunction ^(x, t) obeys the Schrodinger equation 



where m is the mass of the Bose particles and V(x — x') is their interaction. In addition, 
n(x, t) is the particle density which will be determined by an additional mean-field theory 
that, for a supersolid, shows the instability associated with the formation of a density wave. 
Hence the (thermal) average (n(x,t)) is periodic in space and independent of time. As 
a result the condensate wavefunction is, if we neglect density fluctuations, also periodic 
and we have indeed both diagonal as well as off-diagonal long-range order. Fluctuations 
in the density, however, induce variations in the phase of the wavefunction and therefore 
in the superfluid velocity. Because these density fluctuations can be caused by both lattice 
vibrations and oscillations in the net defect density we conclude that both mechanisms 
can lead to superfluid motion. Together with the existence of a conservation law for the 
net number of defects, this explains from a more microscopic view why we found two third- 
sound modes and two modes with an oscillatory net defect density in the case of a supersolid 
adsorbed onto a substrate. 

For the superhexatic phase we have shown that the hexatic long-range order leads to 
an additional (as compared to the superfluid) hydrodynamic degree of freedom that affects 
only the transverse modes and is therefore at long wavelengths decoupled from the superfluid 
momentum density. This can also be understood from the above picture, since variations in 
the orientational order parameter do not lead to density fluctuations in first instance. As 
a result we find on a substrate only two third-sound modes and thus at the hydrodynamic 
level of description nothing to distinguish the superhexatic from the superfluid. Although 
this is in agreement with the experiments of Chen and Mochel, who indeed only observe 
one third-sound branch below the second critical temperature T c , it is unfortunate for the 




(69) 
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purpose of suggesting a possible identification of the superhexatic phase. On the basis of 
our results we can, however, conclude that a more microscopic probe is needed if one wants 
to detect the orientational order present in a superhexatic helium film. In our opinion this 
appears to be an important, but also difficult experimental challenge. 

Finally, we would also like to point out the possible relevance of our results to the recent 
experiments with bulk solid 4 He.i In these experiments Lengua and Goodkind observe at 
sufficiently high frequencies an additional (resonant) attenuation and velocity change of 
sound. Moreover, they notice that their data can be explained by a simple model of two 
coupled wave equations which turns out to be identical to the longitudonal part of our solid 



hydrodynamics derived in Sec. |T|. Because our two-dimensional hydrodynamics should be 
able to describe the propagation of sound perpendicular to the c-axis of hep 4 He, this confirms 
the conjecture of Lengua and Goodkind that the collective mode observed is associated with 
the motion of defects. For a more detailed discussion of the coupling between sound and 
the defects one should of course consider the fully three-dimensional situation and include 
the anisotropy of the hep crystal. Work in this direction is in progress. 
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FIGURES 

FIG. 1. A negative disclination in a hexagonal lattice. It is formed by removing a 60° wedge 
from the lattice, and then distorting the lattice so that the open edges meet. A positive disclination 
(not shown) is formed by the insertion of a 60° wedge of material into the lattice. 

FIG. 2. A Volterra construction for a pair of dislocations: (a) a section of the two-dimensional 
lattice bounded by the loop C is removed and (b) the edges of the loop are pulled together to form 
the dashed line C. The endpoints V of this line correspond to the positions of the dislocations. 

FIG. 3. Dispersion curves for (a) the transverse and (b) the longitudinal modes in the hexatic 
phase. In both cases v = 1/2. 

FIG. 4. A Volterra construction for an interstitial in a hexagonal lattice: (a) three half-infinite 
lines of lattice points are removed (as indicated by the arrows) and (b) the lattice is distorted by 
drawing together the edges about the removed lines in such a manner that the crystal symmetry 
at large distances from the interstitial is restored. The removal of the three half-infinite lines is 
simply the Volterra construction of three dislocations. Thus an interstitial can be viewed as being 
made from three dislocations. A similar construction can be made for vacancies. 
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